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The form and the velocity of solitary, or indefinitely long, waves 
in a single liquid have been examined experimentally by Scott 
Russell and mathematically by Boussinesq and Rayleigh. The 
much wider problem of the possible aperiodic wave forms at the 
common boundary of two superposed liquids does not seem to have 
received similar treatment. Those who have treated the subject 
of waves of finite height at the surface of separation of two liquids 
have dealt rather with the case of periodic waves, for which a 
different method is suitable. (Priestly, Camb. Phil. Soc. Proc., 1910 ; 
Lamb, ibid., 1922; Kolchine, Math. Ann., 1927-8.) 

The discussion here given follows the method used by J. H. 
Michell in unpublished work. 

The motion is supposed two-dimensional, and will be treated as 
steady by choice of an origin of coordinates moving at the rate of 
the wave-form. The axis of x is taken horizontal and the axis of y 
directed upwards. The independent variables are changed from 
x, y to x, # where y is the stream function for the motion. This 
simplifies the treatment of the conditions over the boundaries, the 
coordinate y being constant over each of them. The dependent 
variable to be found in terms of x and y is now y, for which, there- 
fore, a differential equation must be found. When y is found the 
form of a boundary is given in Cartesian coordinates by ascribing 
the corresponding constant value to yw. 

In carrying out the process of approximation we take as the 
general mathematical characteristic of the long-wave motion that 
the variation of a quantity specifying it (in particular, the gradient 
of the wave form), in a distance equal to the depth of either liquid, 
is a small fraction of the quantity itself. Thus, if we take the unit 
of length as of the order of magnitude of the depth of either liquid, 
the second derivative d?y/dx? is to be a small fraction of dy/dx, and 
so for higher derivatives. The assuniption is to include the small- 
ness of dy/dx itself. The general discussion terminates in the 
expression of the gradient dy/dx of the wave form in terms of y. 
I have considered the conditions under which the gradient takes 
the factor form appropriate to either a crested or an inverted 
(trough) wave form. The expression of x in terms of y in general 
involves elliptic integrals of the third kind. Where the undisturb- 
ed depth of the lower liquid is small we may find an approximate 
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equation involving an elliptic integral of the first kind only to deter- 
mine the form of the symmetric wave. I have dealt, finally, with 
aļfcase of asymmetric wave (bore) where the gradient-equation for 
the form can be integrated without further approximation. 


The Differential Equation for y. 
In terms of independent variables x, y, the corresponding com- 
ponents of velocity are given by 
a ay / ay 2 


v= ay/a Me 
and the vorticity by 


w = ay jax? + ay /ay® 


When the independent variables x, y, are introduced we have 
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Whence, substituting for u and v, 
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-that is, 
mee) = 20129)" 5 Oy ey 2 ay)*) 22y 
22) ax? 34) “an ay a. 2) lan 


‘Therefore for irrotational motion, where w=0, we have 


BON 2B Sex HER- 


To investigate a type of irrotational waves we must now find an 
approximate solution of this equation which will satisfy also the 
boundary conditions of the problem. 

J. H. Michell has used this process as an alternative method of 
determining the well-known results for the infinitesimal and solitary 
long waves at the free upper surface of a liquid. The method 
applies equally well to problems on superposed liquids, and I have 
used it to find the equation to the form of the wave of finite height 
and wave length as far as the terms of the sixth order in the wave 
height. 

The question to be considered here, however, is the form of the 
long wave at the boundary between two liquids in relative motion, 
the whole being confined between parallel planes at a distance h 
apart. 

Bist y=0, y=h be the fixed horizontal planes between which the 
liquidshe. Lety=0 at y=0, y=a at y=h and y=c at the interface 
of the liquids. Finally, let p, p’ be the densities of the lower and 
upper liquids and U, V their respective “undisturbed” velocities. 

At the first step in the approximate solution of the differential 
equation (1) for y, we neglect the first two terms as of the second 
order and the equation then reduces to 


2 
EES a S: (2) 
ay? 
On integration, this gives, for the lower liquid, 
AV 70h. aE spe eeee (3) 


where 7 is a function of x. (There is no term independent of y 
since y=0 when y=0.) 

Substituting the value of y given by (3) in the second order 
terms of (1), and integrating again, we find 


2 ama 
y + aye (ne 274 2922] )=m, a re (4) 


ax? 


and putting yyy in the second order terms of (4), we obtain 


ay 2 5 
y } Z = ayes )= nye eeesesoesssossesesoes (5) 
Using the result 
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and its consequence 
=? 44 {522% 9929), 
Wy Vout 3 aT \ 
we can write (5) in the form 


2 2 ay 
27) _ oy? J }=y 2 sitaseeeeneonnens (6): 


an ay 2y 
y Hy ax? OX ay 


For the upper liquid, when we integrate the equation 92y/ay2 =0- 
we get 


since y=h when p =a. 
Following the same steps as in the case of the lower liquid we 
get the equation 
2 2 
j ~4{ 92 2, (22) = - eee 8 
yy BES — Hy -m (B2) b= (y - a) 29 (8) 
Since the pressure must be continuous across the interface, we 
deduce from Bernouilli’s pressure equation the result 
pg? =p == (A Sapp hy sscnecseaecoees (9) 


for points on the interface, where q and g’ are the velocities in the 
lower and upper liquids respectively at the point considered, and 
A is some constant. 


But pets 2 ; 
q Ox 


= ay 
aw 


and y=c, at the interface, so from (6) we find that at the interface 


ay _ H 2232y ay ih 
t — =y-s — — 2 ee 
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aea 9 ay ay ') 
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and in a similar way we find 
are 72 2 
qv= Gantt ae 3 =% az) L he ae (11) 
Hence (9) becomes 
pe? p(c-a)” q 2fpc? _p'(c—a)*\ ay 


ye Q-A | Rly O-R) Sax? 


and therefore 
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2 Ui Lae) ay\2 
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and this is the differential equation for the form of the interface. 
We may write it 


pe? _p'(c—a)? k 2- p'(c - aye (2)° 


pe? p'(c —a)?)\ (dy\2 
-a ~ (y-hyP (2) = (A - 3g) =p) oee (13), 

that is, 

4 cpe?_ p'(c—a)*) dy)? _ E _ oy) 3007, 3 (e—a)? 

ali y—h }(2) ]=304 2g) (p DE ¢ Gone 
eine (14): 

Integrating this we obtain 

pe?_pl(o—a)*\ (dy\*_ ay pray. 4 Bae , 3p! ea)? 

5 = G) ee) Pane TY 
oe (15): 


where D is a constant of integration. 
Thus 


(y= 4 D+By - 3g(e—p)y? p y(y—h) +3pc?(y —A) - 3p' (c —a)?y 
dx pe?(y —h)—p'(c — a)*y 


where 3A (ep—p')=B. 

This is the expression found by J. H. Michell for the gradient. 
We now assume that this expression will factorize in such a manner 
as to give the desired wave form, and then consider the further 
conditions which will make such a form possible. That is, we 
suppose 


dy\* _ —38(o—p')(y ~*)(y —ki)(y — Aa) 
(37) — oe saty a EE i (17) 


This makes dy/dx=0 and d?y/dx?=0 when y=k; and dyjdx=0 
when y= k, and when y= ks. 

Thus with this form the condition that the surface may be 
horizontal when y=4, is satisfied. 

Now for (16) to be equivalent to (17) we must have, by equating 
coefficients of y, 


2 
k?kika = i eee ae (18), 
3g(p—p') 4 kkikatk?(kıt+ka) | = —Dh+3 4pc- p'(c—a)?%, ..(19) 
—3g(p—p') d R + 2k (kı the) thyky =D- Bh,  ..cccecesscccsseeeees (20) 
LOR R FA ERa | B 3gA( op), eeren ea (21) 
and from these equations (18)-(21) we deduce :— 


pc® p'(c—a)? 
a ae — ony eee 
hit ham pay p) -Age * 
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pc*h 2 
= o A? sosucerosecocsesecosooes 23 
and Maham aap) p 
so that /y, ką are the roots of the equation 
pe? p (c — a)? pcrh 
2 Íh ak ea EPLE ON OA 
a pepo! Re (24) 
We therefore have 
2 , 2 
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We may write 


A (c— a} 
seit te Zya 
RB and ay 
since U is the velocity at infinity of the undisturbed lower liquid 
of{depth k, and V is the velocity at infinity of the undisturbed 
upper liquid of depth (A ~ k). 
If we also write p’=Ap and V2= pU? we have :— 


B U2 af fy Ua ae 
oo ys ae AJ {h+ gia)! (gil=ay: 


A aoa Eoo e Ee i 


Necessary Conditions for such a Wave. 
We have put the equation for the gradient into the form 


(y 38 (1-A) (V = RPO -k(O ~ko) | 
a) ~ UR — 1 kA- k) by] 


Now the denominator may be written k? (h - y) +a (4 -)2y and 
y is less than / at all points on the interface. Therefore the denomi- 
nator is always positive. Hence, assuming \<1 (i.e., p’<p), we 
must have y- kı and y - kg of the same sign, to make dy/dx real. 
But y lies between & and either k, or ky, since k, ky, and k, are 
the turning values of y. Therefore, either 
(i) Rey<kh—kg, 
or (li) k>y>k;>kų 
These alternatives represent 
(i) a crested wave, 
or (ii) an inverted wave. 
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There is no wave for a value of & between k, and hg. 

Thus, for values of & between 0 and Å, there is a crested wave, 
and for values of & between ką and + there is an inverted wave. 

Lamb has treated the infinitesimal wave at the interface between 
two liquids (see Lamb’s Hydrodynamics, Arts. 231-234), and if 
in Lamb’s result we make the wave length tend to infinity, we 
find, as we should expect, that the two heights at which infinitesimal 
long waves are possible are k, and ky. 

Now, since k, will be the height of the crest when a crested wave 
exists and +, will be the depth of the lower liquid at the trough in 
the case of an inverted wave, it will be necessary for k, and k 
to be real if there is to be a wave form at all. Therefore, referring 
to the equation (24), we deduce the condition 


UXl—-Aw)2  4hU? 
{a + a gd -A 


Approximation-Method for High Waves. 
If we take & very small we find approximately 


ey c Lia 2k y= Pay — Fa) 
(a y 
= C(y — 2h)(y - hy)(y - ka), 


for values of y near the crest, where C is a known constant. This 
makes 
£ dy 
VA a DA =/ . 
z V(y = 2h) y — Ra) (Y - ka) 


Hence we can find an approximate form for the wave in terms 
of an elliptic integral when the wave is near its greatest{height. 


The Asymmetric Long Wave. 


There is, however, a type of long wave whose form can be deter- 
mined from the differential equation without further approximation. 
This is the wave which we get on putting k=k. Its differential 


equation is 
j- 3g(1 —A)(y - k)? (kı -y 
dx! ~ UJRh- 3 kè- rph- kè by] 


and therefore when y=, or y= k dy/dx=0, and d*y/dx*=0. 

This means that the wave has no crests but rises gradually, 
through an infinite horizontal distance from y=% to y=k,. The 
motion is here of the nature of a “bore.” 

Since 4, and $s are the roots of equation (24), the condition that 
kı should be equal to k; is 


o pea)? tbc 
UF a(p —p) (h—k)®g(p—p’) } Sn a ore (25) 
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that is, 
( U? V2r "i 4U? 5 
Sa ae dT en l r E E A E TAO d 
E l ales) wae 
This gives 
A 2 (U2— ÀV?) 
han aT A O EAE E =; GIOCO D00 (27) 


and the roots of this are always real, since (U?+AV)2> (U2—AV2)2, 
This means that for any given pair of values of the velocities U, V, 
of the currents, there are two possible values of h, the distance 
apart of the horizontal boundaries ; they are given by 


1 7 hi 

h =T — VU 26 AWB, obeancooncenensocnone (28) 

If we regard equation (26) as an equation for V in terms of h and U 
we find 


|\vi= Se e (29) 
When condition (26) is satisfied, we have, from (22) 
re Vo) 
hahah + E E er S.. (30) 
and on substituting for V from equation (29) we deduce 
ae 
= z ser scinsie ddd spinvevaretonge 3 
Ri i oy 
The positive sign with the root in (29) would give 
U2 
oe ae) 


and we consider, therefore, only the negative sign. That is, we take 


ave=U4{1 Ae. 


vulta 
=U K 4 1). 
Now, returning to the equation for (dy/dx)*, these results give 
— 8g(1—A) (y= A)?(Ry - y)? 
dy/dx)? =? ee VE 
Cra U? [M1R h-ki hjk) ty] 
LA (y— k)? (kı —y)? 


Uh [kk 2- 4 kit (h—k)*(h— ky)? + y/h] 


3(y — k)" (kı — 9)? : 
[kk — 4 RPK — (h= k)?(h— Ry)? ¢ yh] 


Solitary Waves. 173 


“Three cases now arise, depending on Ene 

(i). kk,y=(A—k)(h—Rk,), 

(ii) khy> (h—k) (b— Ry), 

(iii) kk, = (h—k)(h—Ay). 
We shall now consider these separately. 

(i) Here we have kk, = (h—k)(h—k,) 
„and, therefore, k-++-A, =h. 
This means that the highest and lowest levels of the wave are 
equidistant from the mean height 4/2 of the liquids. 


In this case 
ye 3(y —k®)(k, —y*) 
k2k,? ‘ 


d 
= O-A) 


where y lies between $ and fy. 
On integrating, if we choose the origin so that the constant of 
integration is zero, we find 


„and therefore 


in? JEA ktk)? 
S 
k+k 
=- 
kak 
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Now changing over to a horizontal axis along the mean level, so 


ithat pe et we find the equation to the wave form is 
a os chee!” A 
1 i= 
that is, y’=a tanh mx, 
where pestis k 
2 
AM N 
sa r 
(ii) Consider now the case kk,> (—k){h— ky). 
Here 
(o 3(y —k)} (ki —y)? 
= [RRP {RRP (hk PA ki)? tyh] 


_ aly — R)? (Ry)? 
~ B? -y 
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3h k?k,? 
j ¢ e ae AA, 
where a ERa Pa and £ ae 
Now, since kk, >(h—k)(h—R,), 
therefore A<k+hy, 
and therefore A-—k<k, and h—k,<k. 
RERAN 
p TAE S 

But O'S BRU 

h 
(h—k)*(h—hy)? 


Ibe 


Phe 
h 
=; where 0<é<1. 
Therefore B?>h, 
ins 
S/o 
Now dx _1} ooh _, therefore 
dy “a (y—k)(ky—y) 


(hy “ns f FD -2 dy +f Po = ay 


= -2y7f?—k artanh E 2 + avk kı arcoth fae = 


which is the equation to the wave form in this second case. 
(iii) If kk, <(A—k)(h—-k,), 


Gia 2 ees 


N C 
WES E y e e 
and p= rae, 
Therefore = VERY and 
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Therefore 
(ky—k)ox =2,/ k, +2 artarth v 
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